Study of the volume and spin collapse in orthoferrite LuFeOa using LDA+U 
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Rare earth (R) orthoferrites RFeOs exhibit large volume transitions associated with a spin collapse. We 
present here ab initio calculations on LuFeOs. We show that taking into account the strong correlation among 
the Fe-3rf electrons is necessary. Indeed, with the LDA+U method in the Projector Augmented Wave (PAW), 
we are able to describe the isostructural phase transition at 50 GPa, as well as a volume discontinuity of 6.0% at 
the transition and the considerable reduction of the magnetic moment on the Fe ions. We further investigate the 
effect of the variation of U and J and find a linear dependence of the transition pressure on these parameters. 
We give an interpretation for the non-intuitive effect of J. This emphasizes the need for a correct determination 
of these parameters especially when the LDA+U is applied to systems (e.g in geophysical investigations) where 
the transition pressure is a priori unknown. 

PACS numbers: 62.50.-p,71.27.+a, 71.15.Mb,77.80.Bh 



I. INTRODUCTION 

The magnetic spin collapse under pressure in transition 
metal oxides has attracted great interest in the past few years, 
not only because of the geophysical implications'-^'^^ and fun- 
damental questions on the origin of the transition, but also 
because experiments and structural calculations have become 
feasible in the pressure regions where the spin transition takes 
place (see for example Refs. HHI^S)- These high pressure 
regions could be reached because of the development of the 
diamond anvil cell (DAC)^ on the experimental side. The de- 
velopment of electronic structure codes allowing for structural 
relaxations facilitates the computational treatment. In partic- 
ular, plane wave methods combined with the Projector Aug- 
mented Wave (PAW) frame wor k^i 10,11,12 ^jjq^^ {q treat atoms 
throughout regions, where the ionic radii, the ionic positions, 
the nominal valence as well as the crystal structure might 
heavily vary. However, the standard treatment using the Den- 
sity Functional Theory (DFT)ii^ in the local spin density 
approach (LSDA) is erroneously cumbered by the so-called 
self-interaction energy and more generally, by the wrong de- 
scription of interactions of electrons inside localized orbitals 
(e.g 3d). 

The limitation of DPT methods has provoked the develop- 
ment of new theoretical methods, such as the LDA+U and 
LDA+DMFT method a'^-"'-'Ai^ (from a combination of the 
DFT in the Local Density Approximation (LDA) and a Hub- 
bard Hamiltonian), the self-interaction corrected LSDA (SIC- 
LSDA)r^ or the hybrid functional methodi^ 

Recently spin and volume collapse isostructural transi- 
tions under pressure have been observed in orthoferrites2ii2^ 
with X-ray-diffraction methods and Mossbauer spectroscopy. 
LuFeOj is an ideal material to test the agreement between ex- 
perimental and computational methods. The reasons for this 
are manifold: First, accurate experimental data exist up to 
pressures of 125 GPaJS Second, the transition is well defined 
contrary to some other orthoferrites (e.g. PrFeOs). Third, this 
compound has a simple magnetic structure,^-' because of the 
complete /-shell of lutetium: correlation effects inside the /- 
shell can thus be neglected and the focus can be put on the 
correct description of the iron atom. Fourth, the distortion 



of the perovskite structure in LuFeOa is strong which clearly 
defines the geometrical structure of the crystal and allows to 
neglect thermal effects on the crystal structure (in other per- 
ovskite materials21 such as SrTiOj and BaTiOa this distortion 
is smaller and leads to a sequence of thermal phase transi- 
tions). This simplifies the theoretical treatment and allows 
straight forward comparison with experimental results. Fifth, 
the stoichiometry of LuFeOa perovskite is well defined in con- 
trast to other materials (e.g. FeO) where it is difficult to obtain 
pure samplesis^is^ 

As a large number of transition-metal compounds are in- 
sulators, the existence of correlated metals raises many the- 
oretical questions which have been tackled recently (see e.g. 
Refs. 27,28i22-) However, orthoferrites are Mott insulators 
with a large gap ^''^~-^° and the rare earth ions retain an impor- 
tant atomic behavior. 

In general, the high spin (HS) to low spin (LS) transition 
in these systems is linked to the considerable volume collapse 
of the transition metal ionAL It leads to the violation of atomic 
first Hund's rule because of the enhanced crystal field and thus 
to the considerable reduction or even complete vanishing of 
the magnetic moment. The success of recent calculation using 
LDA+U and LDA+DMFT to describe volume and moment 
collapse in simple oxides such as MnO^ and CoO?- empha- 
sizes that this transition is clearly linked to the existence of 
strong interactions. 

The electronic structure of LuFeO^ has been studied within 
the PEE formalism by Xing et alM and Iglesias et al.— 
Though, the volume collapse transition was not studied in 
these works. In PB B?^'^\ the AFM magnetic structure is 
found correctly but no band gap^ is found or it is small 
(0.46 eV=). Indeed the R-FeOs orthoferrites are known for 
their large optical gapiai*^ In Ref. [lU atomic relaxations were 
applied whereas in Ref. [s^ the calculations were performed 
in the ideal cubic geometry, which could explain the different 
findings. Our work supports the notion of Ref. [34]: The com- 
bination of a standard treatment of the electronic exchange- 
correlation energy (LDA here, GGA in Ref. [34l) with the cu- 
bic crystal structure results in a metallic state. Recently, Singh 
et al.— have performed calculations using the LDA+U for- 
malism: much larger gaps are obtained. Unfortunately, the 
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FIG. 1: The enthalpy of the HS, the LS and the LS' phase of LuFeOs 
in the pressure range of 0-100 GPa using relaxed (relx) and unre- 
laxed LDA+U calculations. The enthalpy //" (p) obtained from the 
fit of the BM3 EOS was used as a reference. Further details in the 
text. Inset: the energy of the HS and LS phase at the corresponding 
volumes. 



value of Coulomb interaction U is not given. 

In our study we thus focus on the phase transition, and use 
the LDA+U approximation which describes well Mott insula- 
tors such as LuFeOj. We show that the method is indeed able 
to describe the volume collapse associated to the spin transi- 
tion upon pressure. We study how U and / contribute to the 
stabilization of the two phases. 



II. METHODOLOGY 



TABLE I: Summary of the PAW atomic data. The values given 
represent our calculated quantities. The data of Ref.[38| are obtained 
from ultrasoft pseudopotential LDA calculations. 
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points in total, energy cutoff 20 Ha ^ 545 eV), rhombohedral 
iron oxide and hexagonal Lu metal (k-point: 19 x 19 x 13, 
585 points in total, energy cutoff 20 Ha ^ 545 eV). Equi- 
librium properties are compared to values in the literature in 
Tab. U which validates the PAW atomic data. In the case of 
lutetium metal, calculations of Ref. |38i are carried out with 
the ultrasoft-pseudopotential method, which could explain the 
discrepancy with the more precise PAW method, especially 
for a system which contains strongly localized / electrons 
(mainly contained in the augmentation region). 

In order to test the validity of the PAW datasets we per- 
formed LDA calculations on rhombohedral FeO. A four 
atomic rhombohedral {R-3m, space group 166) unit cell 
was chosen in order to accommodate antiferromagnetic 
orderingi^ We use 30 k-points for the k-point sampling and 
a 16 Ha plane wave energy cutoff. Thus aq is converged to 
0.001 A and Bq to 7 GPa. The resulting EOS parameters of 
the AFM phase are given in Tab. U 



B. Calculation setup 

1. Computational scheme 



In this section we gathered the description of the PAW 
datasets, the computational details related to LuFeOs, an anal- 
ysis of the atomic occupation matrices for the HS and LS 
phases in the LDA+U method and a discussion of the ther- 
modynamic potential at pressure, the enthalpy. 



A. PAW atomic datasets 

PAW atomic datasets are generated using ATOMPAWi^^i^ 
For the LDA+U PAW calculations, semi-core states of Lu 
and Fe are treated in the valence. Valence states for Lu, Fe 
and O thus include 5s5p5d4f6s, 3s3p4-s4p3d and 2s2p states, 
respectively. PAW radii are 2.52 a.u., 2.01 a.u. and 1.11 a.u., 
respectively. These values of the radii have been chosen in 
order to avoid the overlap of PAW augmentation regions at 
high pressure. 

The Lu, Fe and O atomic data were tested for the oxygen 
molecule (O2, box size 10 aq - 5.3 A, energy cutoff 16 Ha^ 
545 eV), bcc ferromagnetic iron (k-points: 11x11x11, 56 



Calculations are performed using the ABINIT packager* 
within the PAW— framework. The electron-electron interac- 
tion is treated using the LDA and the LDA+U.i^'i^'^ii^^ 
The unit cell of the ideal perovskite structure is cubic and con- 
tains 5 atoms. In order to establish the GdFeOs orthorhombic 
distortion observed by Rozenberg et al.— a unit cell has to 
be considered with lattice vectors ei + 63, 2e2 and ei - 63, 
where e, are the lattice vectors of the ideal cubic perovskite. 
This unit cell contains 20 atoms and allows for antiferromag- 
netic ordering of the Fe atoms. This unit cell was used for 
our unrelaxed antiferromagnetic calculations while the sym- 
metry elements were reduced to those of space group Pbnm. 
For the unrelaxed structure, we found the antiferromagnetic 
(AFM) configuration to be lower in energy than the ferromag- 
netic (FM) configuration by 1.458 eV in good agreement with 
experiment where the AFM is found for the HS and the LS 
phases,^-^" This AFM symmetry was therefore imposed dur- 
ing the calculation which on the other hand reduced the de- 
grees of freedom. On the Lu positions we find no magnetic 
moment and therefore no magnetic ordering. 

Convergence of LuFe03 was reached on a 6 x 4 x 6 
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Monkhorst sampling grid (36 k-points) and with an energy 
cutoff of 16 Ha («i 435 eV). Energy differences of the HS 
and the LS phase were converged to 10 meV. The pressure is 
converged to 0.08 GPa and enthalpies to 50 meV. Moreover, 
transition pressures are converged to 0.04 GPa. The double 
counting energy in the fully localized limit (atomic limit) 

£dc[«n = U/2N{N-1)-J/2[nHn^ -l) + N^(N^-l)] (1) 

was used throughout this work, where A^"" = tr and - 
+ jyi ^44.45.48 j-^^ parameter U representing the Coulomb 
repulsion of the Fe-3d orbitals was chosen 4.3 eV as it was 
determined by Cococcioni and de Gironcoli- for the oxide 
FeO. A higher value could be used because of the more con- 
tracted orbitals in Fe^+ with respect to Fe^^4^ However, in 
other works-'" a slightly lower value was determined (3.7 eV). 
The values are thus difficult to compare because they depend 
on the orbital basis set. For the sake of simplicity and in order 
to compare with the Ref. the exchange parameters / was 
set to 0.^' Spin-orbit coupling is neglected in these calcula- 
tions. We assume that the orbital magnetic moment - although 
not negligible - is mainly quenched by the crystal field. 



2. Determination of the ground state 

In the ionic limit, the charge of lutetium, oxygen, and iron 
would be 3-I-, 2- and 3+. This implies a formal occupancy 
of 5 electrons for the ^f-orbital sub-shell. In a cubic lattice, 
the ground state would thus consist of the filling of 3 t2g and 
2 eg orbitals (Hund's Rule requires maximal spin polariza- 
tion), while in the LS state electrons fill 5 t2g states. In the 
atomic limit — corresponding to the complete filling of or- 
bitals and without fluctuations — , and thus in LDA+U, the 
symmetry would therefore be broken in the LS stated Calcu- 
lations for both phases have thus been carried out by imposing 
the Pbnm symmetry to the electronic states even for the undis- 
torted structure. The experimentally observed crystal distor- 
tion appears not to be the consequence of the electron distri- 
bution in the LS phase alone, as the distortion is also present 
in the high spin phase. It is probably more due to the geomet- 
ric redistribution of space between the ions of different radii 
as described in section UlI CI 

While for the HS configuration only one electronic ar- 
rangement is possible, the LS can be implemented in several 
ways because in the Pbnm symmetry the degeneracy of all d- 
orbitals is lifted. In order to enhance the convergence of the 
specific spin state the electron-electron interaction potential 
was fixed in the Hamiltonian according to a given occupancy 
matrix (see III for the d correlated subspace during the 
30 first steps of the energy minimization procedure. Then it 
was self-consistently optimized. We determined energies and 
orbital occupancies starting from one hundred different occu- 
pation matrices (for the correlated subspace) at a volume of 
156.92 (^ 90 GPa). These one hundred occupations ma- 
trices correspond to the number of possibilities to have 3 spin 
majority and 2 spin minority states in the ^f-shell. 



The optimization of the electronic density ended in 3 pos- 
sible electronic states which are all tV't^. These 3 electronic 

2g 2g 

States corresponds to different possible coupling of the t2g or- 
bitals. We compared the energies of the 3 electronic configu- 
rations at pressures from to 90 GPa and found that the en- 
thalpies differ by a constant, which is independent of the pres- 
sure. This indicates an identical EOS for the three states in 
the LS spin region (see also FIG.[Tli. Assuming only diagonal 
t2g density matrix as starting point, we find that the configura- 
tion d\^.d^„d^^, is the most stable one.^^^ It is referred to as LS 
in the following. The use of starting density matrix contain- 
ing eg occupations (among the one hundred used) has enabled 
us to study a larger set of solution. Among them, the con- 
figuration LS' has the lowest energy (Els - £^ls' = 0.52 eV 
at a volume of 156.92 A-'). The configuration in the majority 
channel is c/',,ii'.t/J.„. In the minority channel two t2g states are 
occupied and characterized by non-diagonal occupation ma- 
trix elements. It gives a transition pressure which is 4 GPa 
lower than the one of LS. The LS and the LS' structure con- 
verge towards a unique state upon structural relaxation. This 
is reflected in the energy as well as in the occupation matrix 
of the Fe-3ii electrons. The third LS state (LS*) lies between 
LS' and LS, as far as the energy is concerned. In its minority 
channel, only one of the t2g states contributes to non-diagonal 
matrix elements in the occupation matrix. 
Atomic relaxations were performed using the Broyden- 
Fletcher-Goldfarb-Shanno minimization,^ Enthalpies were 
converged to a precision of at least 0.2 meV, energies to 
10"^ meV, forces to 0.01 eV/A. 

3. The thermodynamic potential at pressure 

The 3rd order Birch-Murnaghan (BM3)^i^ equation of 
state 

£(x) = £o + ^-^[B',{x^^' - If + - If ■ (6 - 4x^/^)1 
p{x) = ^(x^/^ - x^/^) . [1 + \{B', - 4)(x2/3 - 1)] 
x= ^ 

(2) 

was fitted to the energy-volume data, where £0 is the 
ground state energy, Vo the ground state volume, Bq the bulk 
modulus and its first volume derivative. 

In density functional theory the correct density is the one 
which corresponds to the lowest energy. Two different phases 
can correspond to local minima of the energy surface E{V), 
where E is the total energy and V the volume. 

However, the free parameter in experiment is not the vol- 
ume V but the pressure p applied to the sample. The cor- 
responding thermodynamic potential is the enthalpy H{p) - 
E(V{p)) + p V{p). The enthalpy can be computed if the pres- 
sure p is determined from the electronic density using the 
Hellman Feynman theorem as it is done in ABINIT. The pres- 
sure at which the transition occurs is given by the equality 
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FIG. 2: Total density of states of from LDA calculations (for expla- 
nation of symbols see caption FIG.[3]l. |(a)| The relaxed LDA structure 
at a pressure of -3 GPa which corresponds to a volume of 197.458 . 
The gap vanishes and the local magnetic iron moments 1.1 /ig. |(b)| 
The relaxed LDA at a pressure of -12.3 GPa which corresponds to 
a volume of 224.882 A^^. The gap is 0.7 eV and the local magnetic 
moments on Fe are 3.56 /ib- 
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FIG. 3: Total density of states of LuFeOs (black) and projected den- 
sities on Lu, Fe and O from LDA-l-U calculations. For Fe the spin 
density is shown, positive values for the minority channel (f ), nega- 
tive values for the majority channel (J,). Lu-/ electrons remain local- 
ized and hardly hybridize. For better visibility these peaks were cut. 
|(a)| The relaxed HS phase at a pressure of 3.2 GPa which corresponds 
to a volume of 210.54 A^. |(b)| The relaxed LS phase at a pressure of 
48.1 GPa which corresponds to a volume of 170.21 A'. 



of the enthalpies. The transition pressure can be equivalently 
given by the well known construction of common tangents on 
E(V) curves. 

In practice, the discrete data set {p,. Hi] of various phases 
can be compared to the continuous enthalpy //, (/?) of an ad- 
ditional phase, which can be chosen to be the enthalpy of the 
fitted E{V) curve with the BM3 expression. The enthalpy dif- 
ference AH^ - - Hr (where s e {HS, LS}) allows to study 
the effects that lead to phase transitions and leave away fea- 
tures which are common to all phases considered (FIG.[T]). In 
particular, the transition pressure is given by the intersection 
of two curves. 



III. RESULTS 

This section presents the description of the volume-collapse 
transition in LDA-i-U. 



A. Electronic properties and densities of states 

We first have performed LDA calculations for the ideal cu- 
bic configuration and have found the resulting DOS to be 
metallic and in particular (i-states are at the Fermi level. Even 
starting from converged LDA-i-U ionic structures and the wave 
functions of the HS and the LS phase, respectively, the LDA 
fails: At small pressures (large volume) it gives a metallic 
DOS with a small magnetic moment of (FIG. |2(a)[ ). 

Only at an extremely large volume (highly negative pressure, 
FIG. |2(b)[ ) a gap opens (0.7 eV) and a considerable mag- 
netic moment results (3.56 hb)- This electronic configuration 
though is not stable at smaller volume (in particular not at the 
volume corresponding to ambient conditions). This empha- 
sizes the incorrect description of correlations. Indeed, rare 
earth orthoferrites are know to be large gap insulators /aii^Si^ 
It is not possible to stabilize two different phases at positive 
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pressure within the LDA. 

As described above (section IIIBb and emphasized before 
(e.g Ref. H), the LDA+U method, which takes into account 
the strong correlations in the atomic limit, is able to describe 
these two phases. The density of states (DOS) of the HS and 
LS phase - at a volume belonging to their individual stability 
range - are shown on FIG. |3(a)| and |3(b)| The main effect of 
LDA+U is to stabilize the HS phase at positive pressures and 
to increase the gap inside the tZ-orbitals. In the LS and the 
HS phases, the Fe-3d orbitals in the valence band are strongly 
hybridized with the oxygen p orbitals. The Lu-f retain their 
atomic-like character while the Lu-d - though unoccupied - 
hybridize with O-p orbitals. The DOS scarcely changes with 
pressure, even at low pressures (0 GPa) where the LS phase 
is unstable. The HS phase (FIG. |3(b)[ ) is characterized by 
atomic-like peaks for all elements, especially Lu and Fe. 

At GPa {U = 4.3 eV and 7 = 0, no relaxation) the HS 
phase is insulating with a gap of L75 eV, and the LS with a 
gap of L16 eV. At 50 GPa the gaps are LIO eV and 0.93 eV, 
respectively. In fact, the gap size shows a perfect linear behav- 
ior in the range between and 200 GPa. Relaxation increases 
the band gap at GPa to 2.17 eV (HS) and 2.03 eV (LS). In 
Ref. the value of the gap of the HS phase obtained in 
LDA-hU is ^ 8 eV. The comparison with our results remains 
difficult because the parameters U and / were not given there. 

One of the major impacts of the parameters U and J is their 
influence on the gap size. While U always increases the gap, J 
decreases it. The effect on the HS phase is half as strong as on 
the LS, where the gap increases by 0.45 eV when U increases 
by 1 eV and decreases by 0.85 eV when / increases by 1 eV. 

The local spin moment on the Fe atoms is found to be 
4.05 jUs in the HS and 0.99 //b in the LS structure at GPa, 
which decreases to 3.83 fiB and 0.91 jUb at 100 GPa, respec- 
tively (100 GPa con-esponds to V =164.79 A in the HS and 
y =155.00 A in the LS phase). 

The discrepancy between the ionic ideal value of the spin 
moment and the actual value can be explained as follows: The 
hybridization of Fe-d electrons with oxygen-/? electrons gen- 
erates small occupancies on orbitals which are empty in the 
ionic picture, eg orbitals are more hybridized with oxygen-/? 
states: they are thus more filled. In the HS phase, 2 of the 
hybridized eg contribute to the minority spin and thus to the 
considerable reduction of the magnetic moment with respect 
to the ionic value (i.e. 5 yUg). In the LS phase, the hybridized 
eg orbitals appear in both spin channels and cancel each other 
Only an almost empty t2g orbital reduces the local magnetic 
moment, which therefore remains close to the ionic value. 

As will be detailed in the next section, at ambient pressure 
we describe an insulator-insulator transition for LuFeOs. Re- 
cent LDA+U calculations on MnO give a similar conclusions, 
although the LS state obtained in this work seems to be less 
intuitive.^ However we emphasize here the limitation of our 
work which does not contains fluctuations. These fluctua- 
tions could easily make the system metallic as has been shown 
recently in MnO.-- Moreover, recent experimental works on 
NdFeOs^'' and BiFeOj^^^ show that the transition is closer to 
an insulator-semiconductor or an insulator-metal transitions. 
Further optical experiment or LDA+DMFT calculations on 



LuFeOs could clarify this issue. 



B. The equation of state and the transition pressure 

As emphasized before, the LDA is not able to describe the 
occurrence of two phases. Nevertheless, a continuous and lin- 
ear decrease of the local magnetic moment from 1.05 /ib to 
0.55 //B is observed between and 100 GPa (relaxed struc- 
ture). 

LDA+U calculations for the two phases were first per- 
formed with the frozen ideal cubic configuration on the 20 
atoms super-cell. The symmetry was reduced to the 8 ele- 
ments of the corresponding Pbnm space group in order to al- 
low lifting of the degeneracy of the Vs-d states. The energy 
and enthalpy versus volume curves are reproduced on Fig [T] 
One can see a phase transition between the LS and the HS 
phases. It occurs at 22 GPa. 

This pressure shifts to 5 1 GPa after structural relaxations, 
which compares well with the experimentally observed transi- 
tion pressure of 50 GPa (FIG.© -- Note, that the HS-LS tran- 
sition was not found in earlier works and the Fe magnetic mo- 
ment was underestimated (3.6 ij-b)^- The calculated collapse 
of the ground state volume at the transition amounts to 6.0%, 
compared to the measured 5.5%. Whereas in LDA the volume 
at zero pressure is underestimated by 12%, the introduction of 
interactions localizes the electrons, and thus raises the vol- 
ume: the volumes at zero pressure are then well reproduced 
(within 2.5%, see Tab. for both the HS and LS phases. 
Bulk moduli are in the range of experiment but with an er- 
ror of about 20%. As a consequence, the difference between 
theoretical volume and experimental one increases to 6% at 
150 GPa (see FigHJi. Besides, the discontinuity of the bulk 
modulus at the transition is qualitatively described (Tab. [Till: 
3 1 GPa (this work) compared to 72 GPa in experiments^ (see 
also FIG.©. 

Positional parameters are in good agreement with the ex- 
perimental findings, moreover, we find no significant discon- 
tinuity of the positional parameters at the transition pressure 
as observed experimentally (FIG.|4]inset). 

We investigated the effects of the Hubbard part of the 
Hamiltonian on the transition pressure. We therefore com- 
puted energies (non-self-consistently) in the LDA frame from 
the densities obtained in LDA+U (£'lda[«lda+u]) and ob- 
served that the HS-LS transition is maintained but shifted to 
much lower pressures. This shows that the LDA+U potential 
creates two different electronic densities which give different 
LDA energies. The £u term in the energy additionally de- 
termines the difference of energy between the two configura- 
tions, LS and HS (see also the discussion in section HIl Di ). 

Finally, we mention that we have neglected the entropic 
contribution in these calculations. In the iso-structural transi- 
tion in cerium, the entropy appears to be essential to describe 
the transitioniS In LuFeOa we could expect the variation of 
entropy to be rather weak also because of the broken symme- 
try of d states. However, we could miss fluctuations between 
configurations in the LS case. Furthermore, theoretical stud- 
ies beyond LDA+U (e.g LDA+DMFT) as well as experimen- 
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TABLE II: Parameters of the fitted S"* order Birch Murnaghan EOS 
of the relaxed HS and LS phase of LuFeOs LDA+U ( (/ = 4.3 eV and 
J = 0) . -Eo.HS ~ Eoi^s =3.25 eV. The gap at P =0 GPa is also given. 





Vo [A3] 


So 


So [GPa] 


Gap [eV] 


LDA 




195.06 


4.19 


236 





LDA+U 


HS 


213.59 


3.80 


214 


2.17 


LDA+U 


LS 


197.41 
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2.03 


Expis^ 


HS 


218.40 
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ExpS 


LS 


199.40 
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FIG. 4: The volume of LuFeOs at pressures between and 120 GPa. 
Bullets (•): Calculated values using LDA+U with parameters 
U = 4.3 and 7 = eV and performing structural relaxations. Dashed 
lines: Calculated values using the same parameters U and J with- 
out structural relaxation. The calculated HS-LS transition occurs at 
21 GPa. Boxes (□): Experimental data from Ref.|22. Inset: Reduced 
atomic coordinates of Lu, bullets (•) calculated, boxes (□) experi- 
mental data from Ref . The symbols correspond to the X and Z 
coordinates, respectively. 



tal studies of the transition as a function of temperature could 
help to understand these issues. 



TABLE III: Atomic Wyckotf positions of LuFeOs in the HS and the 
LS phase (space group Pbnm, no. 62). The lattice parameters of 
the HS at GPa are 5.156 A, 7.491 A and 5.530 A compared to 
the experimental values of Marezio et al.— of 5.213 A, 7.565 A and 
5.547 A (maximal deviation for a by -1.1%, which shows that the 
chosen U parameter cannot fully correct the known overbinding be- 
haviour of LDA). The lattice parameters of the LS at 50 GPa are 
4.617 A, 6.896 A and 5.323 A. Compared to the data pubhshed by 
Marezio et al.— for the HS phase the average positional deviation 
is 0.013 A, and the biggest deviation is for Lu by 0.03 A. The unit 
cell volume is underestimated by 2.4% in the calculations^ while 
the experimental ratio of the unit cell vectors a:b:c of 0.689:0.733:1 
is reproduced by 0.688:0.738: 1 in the calculation. 



Atom 


Wyckof position x 


y 


z 




HS at GPa, this work 


Lu 


4c 


0.0233 


0. 


0.0758 


Fe 


4a 


0. 


0. 


0. 


1 


4c 


0.3783 


1/2 


0.0439 


02 


8d 


0.8126 


0.3112 


0.1930 




HS experimental, Ref. 64 


Lu 


4c 


0.01997 


0. 


0.07149 


Fe 


4a 


0. 


0. 


0. 


O 1 


4c 


0.38010 


1/2 


0.04610 


02 


8d 


0.81070 


0.31210 


0.19290 




LS at 50 GPa, this work 


Lu 


4c 


0.0371 


0. 


0.0854 


Fe 


4a 


0. 


0. 


0. 


1 


4c 


0.3941 


1/2 


0.0300 


02 


8d 


0.8267 


0.2987 


0.1917 



The tilting angles are only approximate, because a symmetry 
conserving distortion is superimposed to the tilting of the 
octahedra and in the distorted octahedra the determination 
of the tilting angle is not unique. While in the HS phase 
the octahedra are elongated by 1 % in the LS phase they are 
squeezed by 3%. The relaxed atomic positions of the HS and 
the LS phase can be found in Tab. |III1 The more important 
distortion of the octahedra in the LS spin phase could be 
attributed to the existence of a Jahn Teller effect in this phase. 



C. The ground state structure 

Perovskites (ABX3) are known to consist of mostly rigid 
BXe octahedra, while the A cations are placed in the intersti- 
tial space between the octahedra. Tilting of the mostly rigid 
octahedra allows to decrease space assigned to the A cations 
and the unit cell volume. The tilting angle is sensitive to the 
ratio of the B-X bonding length and the ionic radius of Ai^Si^ 
According to Glazer— and Woodward— perovskites can be 
classified according to three tilting angles and the phase (+/-) 
of successive octahedra along the tilting axis (see Ref. [62 
p. 3386 for a sketch of the tilting system and Ref. |63| p. 34 
and p. 36 for further details on this concept). LuFeOa belongs 
to the space group Pbnm and the tilting system is a^b'b^ . 
We find tilting angles of a ^ b ^ 13° in the HS phase at 
3 GPa and fl ^ 17 ° and ^ 11 ° in the LS phase at 87 GPa. 



D. Dependence of the transition pressure upon U and J 

A number of calculations were performed with varying 
Hubbard parameter U and exchange parameter J. The aim 
was to determine the dependence and the sensitivity of the 
transition pressure on these parameters. We chose the param- 
eters in the reasonable range of [in eV] U — {2, 4, 6, 8) and 
J - {0, 0.7, 1.7}. The transition pressures were calculated us- 
ing the frozen ionic configuration in the cubic structure: we 
were only interested in qualitative trends. The results are vi- 
sualized in FIG. |5] The transition pressure p„ in the region 
explored can be recast as 

p„^A + BU + CJ (3) 
with A = -16.2 GPa, B = 8.75 GPa/eV and C = 
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-10.7 GPa/eV. This expression was obtained from a least 
squares fit to the calculated transition pressures and its reli- 
ability is +1 GPa in the range studied. The resulting coef- 
ficients B and C are comparable with opposite sign, which 
emphasizes that the main parameter of the calculation is the 
difference U - J. It shows that the formulation of Dudarev 
et al.— is a good approximation in the present case. In these 
unrelaxed calculations the increase of U from 4 eV to 8 eV 
raises the amount of the volume collapse insignificantly from 
7.53% to 7.61%. 

The positive dependence of on U can be explained as 
following: Consider the simplified atomic Hamiltonian 



E,,^uY,n]n)+{U-J)J]n>'J 



I, J 

El] — Eee ~ E^Q 



a-,i>j 



(4) 



where the density matrix is diagonal («^ = "^^'7^ ^<ic 
is chosen in the atomic limit (Eq. [1]). E\j can be recast as 
^(N - YjiA^D'y As emphasized before,^ this quan- 
tity cancels when occupation numbers n"' are integers and is 
positive elsewhere (see FlG.l of Ref. l66l). At this point, we 
emphasize that the effect of J is thus not obvious: An increase 
of / stabilizes the high spin state both in E^g and in E^c, so 
that a clear effect on ^tot cannot be simply anticipated as em- 
phasized before for MnO.— 

When going from HS to the LS hybridization effects are 
enhanced. This is partly due to the decrease of the volume 
but also due to the different spin configurations of the HS and 
the LS phases as they were presented in section 1111 Al where 
the reduction of the magnetic moment of the HS phase was 
discussed. As a consequence, c/-orbitals are increasingly hy- 
bridized with p-orbitals. It implies that Bloch states have a 
mixed O p-Fe d character: The d states which should be 
empty in the pure ionic picture (see insets in |3(a)| and |3(b)| i 
are more filled in the LS phase than in the HS spin phase 
(the occupancy is still lower than 0.5). As a consequence: 
E\f > E^^ (see FIG. 1 of Ref. HI). Thus, if U is increased 
or J is decreased, the LS phase is destabilized with respect 
to the HS, and thus the transition pressure increases. This is 
what is observed in FIG. |5] This effect shows that U and J 
cannot be taken as parameters in the calculation, because the 
double counting expression - though approximate - should 
correct the LDA energy. This emphasizes the need for a cor- 
rect determination of U and J. Additionally the effect of J is 
slightly more important than the effect of U. This is due to the 
difference between expression |4] for Eu and the rotationally 
invariant expression that we useii^^ 

The increase of J induces a charge redistribution, which 
further lowers the relative energy of the LS phase (see also 
FIG.|5]insets), i.e. when J increases in the LS phase, electrons 
are transferred from the fully occupied states to the weakly 
occupied states. This transfer can be taken into account in the 
calculation of the LDA+U contribution in Eq.|4] We calculate 
this contribution for a fixed occupation at different values of 
/ and compare it with the energy shift for different values of 
J and relaxed atomic occupancies. We find that in the second 
case the energy change with J is three times bigger than in 
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FIG. 5: The dependence of the transition pressure upon the Coulomb 
repulsion U and the exchange parameter J for the unrelaxed struc- 
tures. An increase of U always increases p„ while an increase of J 
monotonically decreases p^, . Insets: the dependence of the ground 
state energy on U and 7 at a fixed volume, a.) at 174.41 and b.) 
at 183.02 A\ 



the first case. On the other hand, the energy change due to 
electron redistribution is negligible in the case of the increase 
of U. The total number of localized electrons varies when 
U and J vary. This effect is an order of magnitude smaller 
than the two other contributions (i.e. P' the change of Eu 
througth an increase of [U - J) with fixed occupancies and 
2"'' the change of Eu through a charge redistribution with a 
fixed total number of electrons). 

As mentioned before, the bulk modulus is underestimated 
by about 20%, which increases the deviation in volume from 
2.5% to 6%, when pressure increases from to 150 GPa (see 
FIG. Ell. Our analysis of B(){U - J) shows a negative corre- 
lation between U — J and Z?o (with a slope of -30 GPa/eV), 
indicating that the value of U - J was slightly overestimated 
in our structural calculations. In LuFe03 the electron-electron 
interaction is more efficiently screened than in FeO, where the 
interaction parameter was originally derivedi^ 

It is possible, that metallization is the origin of the hard- 
ening of the material at high pressure. As was observed for 
MnO^ and BiFeOs^S the HS-LS transition can be accompa- 
nied by a insulator-metal transition.^- The experimental obser- 
vation could be better described by DMFT calculations^ The 
increased bonding due to metallicity could make the material 
somewhat harder. 
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IV. CONCLUSIONS 

We have carried out LDA+U calculation within the PAW 
framework on the rare earth perovskite LuFe03. We describe 

13 I 

an iso-structural phase transition from a high spin phase f2g^g 
towards a low spin phase fjg^g ^^^^ ^ volume collapse of 6.0% 
(Exp: 5.5%). Atomic positions, magnetic moments and lattice 
constants are computed and are in good agreement with exper- 
imental data.<22 At high pressure, the disagreement on volume 
is at most 6%. The observed reduction of the local magnetic 
moment on iron is ^ 3 fis- 

We find, that the LDA-'?'^- is not apt to treat LuFeOa. The 
LDA+U calculations presented here are always superior be- 
cause the band gap, the phase transition pressure and the local 
magnetic moments could be determined more correctly. 

As the computation of U is not the goal of this work, we 
check the effect of U and / on the transition. We compare 
the filling of orbitals in the HS and LS phases and propose 
an interpretation for the non intuitive effect of /. The deter- 
mination of the parameter of the Coulombic on-site repulsion 
U and the exchange energy J appears to be essential, because 
the critical pressure for the spin collapse depends linearly on 



them. They enter the expression for the critical pressure with 
opposite signs but the same magnitudes [Eq. [3l as shown us- 
ing the simplified LDA+U formalism of Dudarev et al.— 

These calculations open the way to other complex systems 
such as orthoferrites where correlation effects are important 
both on iron and on the rare earth atom. Concerning LuFeOs, 
some improvement of the understanding of the transition 
could also be brought by experimental studies of the tran- 
sition as a function of temperature. Experiments on optical 
properties as a function of pressure as well as LDA+DMFT 
calculations could further clarify the nature of the transition. 
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